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Abstract 



The radiative corrections to the decay processes of the neutral (CP-even) Higgs 
boson (H) into a longitudinal gauge boson pair, i.e., H — > Z L Z L and H — > W^W£ 
are analyzed in the two-Higgs doublet model by assuming that all of the Higgs 
boson masses are much greater than the W and Z bosons'. These calculations are 
motivated to see if one could see potentially large virtual effects to these decay 
rates due to the charged and CP-odd neutral Higgs boson masses (mg and tjia, 
respectively) which are supposed to be larger than m H . It is pointed out that, 
although the radiative corrections to the decay width T(H — > W^W£) depend 
sensitively in general on mc and m^, there occurs a screening effect, i. e., cancellation 
in leading terms once we set mc = rriA, so that the radiative corrections tend to 
be minimized. It is also pointed out that the decay rate T(H — > ZlZl) is fairly 
insensitive to the other heavier Higgs masses and is possibly a good measuring 
tool of the Higgs mixing angle. The mechanism of these screening phenomena in 
the Higgs decays is explained on the basis of a new screening theorem, which we 
postulate with reference to the custodial symmetry in the Higgs potential. 
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1. Introduction 



It is becoming more and more important to investigate non-decoupling effects of 
heavy particles in the low energy observables, since the available energy of future 
accelerators will increase only little by little. One of the examples of low energy 
manifestations of heavy particles has been given by the top quark, whose mass was 
surmised before its discovery [1,2] by detailed analyses of the electroweak data. 

Now that the existence of the top quark has been confirmed, the next pressing 
experimental task is the discovery of the Higgs boson. (For a review, see Ref. [3].) 
It has been known for long time that the Higgs boson in the standard model is very 
elusive: indirect signatures of the Higgs boson appear in the low energy data on 
the oblique-type radiative corrections at most in logarithmic terms at the one loop 
level. This fact is often referred to as Veltman's screening theorem [4]. Perhaps 
more precise electroweak data would change the situation and the signatures of the 
standard model Higgs boson might be just around the corner. In any case, once the 
existence of the standard model Higgs boson would be established, we should be still 
looking after new unknown heavy particles even further, via radiative corrections in 
order to probe what lies beyond the standard model. 

In our previous publication [5] , we have investigated the possibility that one could 
get hold of signatures of unknown heavy particles through radiative corrections, 
supposing that the standard model Higgs boson H has been discovered. As an 
example we adopted a two-Higgs doublet model, the most conservative extension of 
the standard model. Here the Higgs sectors consist of, besides the neutral H boson, 
another (CP-even) neutral Higgs boson h, (CP-odd) neutral Higgs boson A, and 
the charged Higgs boson G ± . Assuming that their masses m#, m^, rriA and m G 
are all much greater than Mw, we computed the radiative corrections to the decay 
width T(H — > W£W£) as a function of nih, tua, and ma- Here the subscript L 
denotes the longitudinal polarization . The preliminary numerical calculations in 
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Ref. [5] show that the magnitude of radiative corrections depends rather sensitively 
on the choice of m G and m A . It is also suggested implicitly that, if uic = m,A, 
then the radiative corrections are minimized. Note in this connection that the Higgs 
potential in the two-doublet model respects an SU (2) L x SU(2) R symmetry [6], if we 
put rriQ = vtia- The existence of the additional global SU(2)r symmetry would lead 
to isospin symmetry SU(2)v = diag[S77(2)z, x SU(2)r] after spontaneous symmetry 
breaking. This isospin symmetry is often called custodial symmetry in literatures. 

The purpose of the present paper is two-fold. The first one is to extend our 
previous analyses to another decay mode H Z L Z L and to compare its dependences 
on ma, uia, uih and Higgs mixing angles (a and /3) with those in H — > W^W^ . We 
would like to establish the following two peculiar facts as to the radiative corrections 
by numerical calculations on computer. 

(1) Although the radiative corrections to the decay H — > W^W£ is in general 
sensitive to the choice of mc and tha, they tend to be minimized if the mass 
parameters are such that the custodial symmetry is respected, i.e., ma = rn A - 

(2) The radiative corrections to the decay H — > Z L Z L are insensitive to the choice 

of the mass parameters. It is only sensitive to the Higgs mixing angles. 

The second purpose of the present work is to give theoretical explanations to (1) 
and (2) without recourse to numerical methods. We will postulate a new screening 
"theorem" for the Higgs decay vertices, by which we are able to acquire a clear 
grasp of tog--, rn A - and m^-dependences of the decay widths T(H — > W£W£) and 
Y{H^Z L Z L ). 

The screening theorem for the Higgs vertices that we just mentioned is in close 
analogy with the celebrated Veltman's theorem [4], which may be applied, in con- 
trast, to the oblique-type radiative corrections. The statement in (1) referring to 
the custodial symmetry reminds us of the detailed study of the Veltman's theorem 
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in the standard model by Einhorn and Wudka [7]. They argued that the radiative 
corrections to the gauge boson propagators may be classified into several types. On 
the basis of the custodial symmetry in the weak [/(l)y-coupling limit, they claimed 
that the Higgs mass dependence at the L-th loop is at most (m#) L_1 rather than 
(m 2 H ) L . The statement (1) indicates that there is a similar situation in the two-Higgs 
doublet model as well, in radiative corrections to the (non-oblique type) decay pro- 
cess if the custodial symmetry is respected. We will make this fact crystalline in 
the form of a new screening theorem. Although we will not go so far as to give as 
general an argument as Einhorn and Wudka's including all order corrections, the 
similarity lying between Veltman's theorem and our counterpart provides a strong 
evidence in favor of the validity of our version beyond the one-loop level. 

Throughout the present paper we will often refer to and make use of some of 
the formulae given in Ref. [5] with respect to the decay H — > W^W£ . In Ref. 
[5], in passing, the effect of the top quark mass was not taken into consideration. 
Since the top quark mass might produce non- negligible effects to the decay [8], we 
will improve our previous analysis by including quark loops as well. Partial lists of 
earlier works on these decay processes in the standard model and supersymmetric 
models are given in Ref. [9-12]. A brief summary of the present work is found in 
Ref. [13]. 

The present paper is organized as follows. Our Higgs potential is specified in Sec. 
2, thereby explaining the connection between the custodial symmetry and the mass 
degeneracy ma = tua- We will make an extensive use of the equivalence theorem 
[14-17] at the loop level. Sec. 3 is devoted to clarification of an issue that arises 
if one uses the equivalence theorem at higher orders [17]. One loop calculations of 
the decay processes, which will be presented in Sec. 5, are rendered complicated, 
because of the mixing between those of same quantum numbers. The method of 
renormalization of the Higgs mixing angles and the wave function renormalization 
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constant matrices are explained at length in Sec. 4. Our numerical calculation 
is presented in Sec. 6. Our new screening theorem is given in Sec. 7 which, we 
believe, will lay cornerstones for clear understanding of the qualitative features of 
the numerical calculations in Sec. 6. Sec. 8 is devoted to summary and discussions. 
Some details of our calculations are relegated to several Appendices. 



2. The Higgs potential and the custodial symmetry 

Before launching into the details of our calculations, we have to specify our 
Higgs potential of the two Higgs doublets, $i and <2> 2 with Y — 1. The criterion 
of determining the Higgs potential is the natural suppression of the flavor- changing 
neutral current [18] and it is often assumed that the discrete symmetry $ 2 — > —$2 
is respected except for soft terms, namely, 

v($ u $ 2 ) = -i4 1 $! 1 2 -i4 1 $ 2 1 2 -(^ 2 2 $I$ 2 + 

+Ai I $1 | 4 +A 2 I $ 2 | 4 +A 3 I $1 | 2 | $2 r 

+A 4 (Re$l$ 2 ) 2 + A 5 (Im<l>l$2) 2 . (1) 
This potential is general enough to encompass supersymmetric models and the soft 
term (/i 2 2 $i $2 + lA^t^i) plays an important role there. The complex phase of 
/ii2 may also be instrumental to the idea of baryogenesis at the electro-weak scale 
[19]. The existence of this term, however, makes our analysis too complicated to get 
an insight into radiative corrections. We will therefore set = throughout our 
calculations. Perhaps it is worth mentioning as another excuse for setting /i 12 = 
that, since we are interested in the non-decoupling effects caused by strong quartic 
couplings (times v 2 ~ (246GeV) 2 ), the mass scale /lu as opposed to v does not have 
direct relevance to what we are concerned with and may be neglected at the first 
step. 
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From the viewpoint of the custodial symmetry, it is convenient to introduce a 
2x2 matrix notation of the Higgs fields, i.e., <3>j = (iT 2 <3>*, <3>j). In terms of this 
notation, our Higgs potential is expressed as 



-^?tr(*f*x) - ^tr($}$ 2 ) 

2 



V($l,$ 2 ) = 

+ ^ {tr^J^) j' + l -\ 2 {tr(*J* 2 ) 
+iA 3 tr($J$ 1 )tr($ 2 t * 2 ) + ^A 4 |tr($J$ 2 ) + tr($}$i)} 2 
-^A 5 |tr(r 3 *|*i) - tr(r 3 $I$ 2 )} 2 . (2) 

This shows clearly that, without the last term in Eq. (1), the potential would pos- 
sess the global symmetry SU(2) L x SU(2) R , under which <frj undergoes the trans- 
formation $j — > gi^ign £ SU(2)l, Qr G SU{2)r). The isospin symmetry 
SU{2) V = diag[S77(2)£ x S , t/(2) R ], that would survive the spontaneous symmetry 
breaking, is also broken by the A5-interactions exclusively. We will come to this 
point later in connection with the screening theorem. 

The particle content of the scalar sector may be seen by putting 



V V2 



; 'vi + hi + izi 



(3) 



and by shuffling Eq. (1). Here v^s (i=l, 2) are the vacuum expectation values. The 
mass eigenstates are obtained by diagonalizing the quadratic terms in the neutral 
as well as charged sectors via 

' h 
H 

■ 



cos a — sin a 
sin a cos a 



cos {3 — sin j3 
sin (3 cos (3 

cos f3 — sin (3 
sin (3 cos (3 




(4) 



(5) 



(6) 



The mixing angle /3 is simply given by tan/5 — vijv\. As we can see easily, H and h 
are CP-even neutral field, while A a CP-odd neutral one. The charged Higgs field 
is denoted by G ± , and the Nambu-Goldstone bosons are w ± and z. 

The physical parameters of our theory are the masses m H , m h , m G , m^ and the 



vacuum expectation value v — y v\ + v\ together with the mixing angles a and (3. 
The quartic couplings in (1) are expressed in terms of these physical parameters by 

Ai = - 2 1 (m^cos 2 a + m^sin 2 a), (7) 
2v z cos z p 

^2 = n 9 1 o o ( m h a + m H COs2 a ) ) ( 8 ) 

sin 2a . 9 9 . 2m 2 -. . , 



2m 



a 4 = — =5 (10) 

A 5 = ^(m 2 4 -m 2 ; ). (11) 

As wee see in Eq. (11), the deviation from the mass degeneracy m 2 G = m\ between 
charged and CP-odd neutral Higgs bosons thus measures the custodial symmetry 
breaking. 

The quartic couplings are all assumed to be real and there is no source of CP- 
violation in the Higgs potential. The CP-invariance enables us to set up several 
selection rules for triple and quartic Higgs couplings. It is worthwhile mentioning 
that we could assign G-parity in the usual way as in hadron physics in connection 
with the isospin symmetry. The neutral h and H bosons have even G-parity, while 
G, A, w and z odd G-parity. This G-parity is also useful to set up selection rules 
for the Higgs- Goldstone interactions for the case of A 5 = 0. 
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3. The equivalence theorem at the loop level 



The Higgs boson decay into a gauge boson pair is dominated preferentially by 
those into longitudinally polarized ones if m# ^> Mz, Mw- In such a case we may 
take an advantage of the equivalence theorem [14-17]. This theorem states that the 
S-matrix elements associated with longitudinal gauge bosons are approximated by 
those of corresponding Nambu-Goldstone bosons; 



T(Z L ( Pl ), Z L (p n ), W L ( qi ), • • • , W L (q m )) 
= (CL A T(CZ d rT(iz {Pl ), ■ ■ .,i Z(Pn ),M<h), ■ ■ -,iw {qm )) + 0(M w /y/s). 



Here y/s is the typical energy scale characterizing the scattering process and C^ od , 
and C^ od are the so-called modification factor [17] to be attached to each external 
gauge boson line of Z and Ws, respectively. 

The equivalence theorem was first proved on the tree level [14, 15]. Since then, 
the validity of this theorem on the loop level has been examined by several authors 
[16-17]. The point is that the right-hand side of Eq. (12) is unphysical (gauge- 
dependent) matrix elements and that we have to specify how to renormalize the 
external lines of the Nambu-Goldstone bosons. The modification factors C^ od and 
C^ od are thus introduced to match the external line renormalization to the physical 
S-matrix on the left-hand side. In a nutshell, if we work in the on-shell renormaliza- 
tion scheme in the Landau gauge, these modification factors turn out at the one-loop 
level to be [17] 



Here the wave function renormalization constants of the gauge bosons (Z and 
W^) and the Nambu-Goldstone bosons (z and w ± )are denoted by Zz, Zw, Z z and 
Z w , respectively. The origin of Eq. (13) may be explained in the following very 



(12) 




(13) 
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intuitive way. The factor ^ Zz/Z z and \J Zw / Z w compensates the difference in the 
external-line renormalization between Nambu-Goldstone and the gauge bosons. The 
presence of <5Mf and 5Myy in Eq. (13) is due to the fact that gauge boson masses 
appearing in the longitudinal polarization vectors in the full theory are simply of 
kinematical origin irrelevant to renormalization, while they appear, in the Higgs- 
Goldstone systems, as coupling constants that are subject to renormalization. Thus 
the presence of 5M| and SM^ in (13) may be undestood as making up this differ- 
ence. 

Toussaint [20] once calculated SM W and <5Mf at the one-loop level in the two 
Higgs doublet model, assuming that all the Higgs masses are much greater than the 
gauge bosons'. For completeness, we would like to include the effects due to the top 
quark mass since it is of interest to see whether the largeness of the top quark mass 
might be just on the verge of affecting the Higgs decays. After including the top 
quark mass effects, Toussaint 's results on the vector boson mass renormalization are 
modified into 

_ 1 \\, 2 , 2 , 2 



M\ (4tt)V\2 



[m H + m h + m A ) 



2 , m 2 A m 2 H , m 2 A . 2 . n . m 2 A m 2 h , m 2 A 

+ cos 2 (a - 0) ^ 9 ln — r + sm 2 (a - 0) ~ \ m ~4 
m H — m A mjj m L h — m A 

+2JVcm?( _±__ T£ _ ln ^_)}, (14) 

8Mw 1 f 1 /n 2 9 9 9 x m 2 r m 2 A . m 2 r 

~nf = n X9 A o ( 2m G + m H + m l + m A + o 2 ln ~ r 

(4vr)V[2 v G " ft A/ - 

9/ ^ rn^rrii, , m 2 ^ 9/ m^rnl , m 2 ^ 

+ cos 2 (a - /3) 9 G H 9 In — f + sin 2 (a - /?) = G \ ln — £ 



+ 2iV cm? (^ + I- 7E -.„^. (15) 



Here iVc is the number of colors (iVc = 3), fi is the mass scale appearing in the 
.D-dimensional regularization method, and here and hereafter we neglect the bottom 
quark mass effects. 
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To compute the wave function renormalization constants Z z and Z w of the 
Nambu-Goldstone bosons and in particular the quark contributions to them, we 
have summarized in Appendix A and Table 1 the Yukawa couplings of the top and 
bottom quarks to various scalar particles. There are two models of the Yukawa 
coupling, so-called model I [21] and model II [22]. In model I, both top and bottom 
quarks receive their masses from only one of the two Higgs doublets, say, $2- In the 
model II, on the other hand, the top quark mass comes from <£> 2 , and the bottom 
quark mass from $i. As we see from the list of the couplings, there exsits little 
difference between these two models as far as the bottom quark mass is negligibly 
small. Setting m b = 0, we will proceed hereafter without dicriminating the two 
models. 

The wave function renormalization constants Z z , and Z w of the Nambu-Goldstone 
bosons are also easily extracted from the calculations in Ref. [20] together with those 
in Appendix C, namely, 

Z z = l-^^^(m 2 H + m 2 h + m 2 A ) 

+ cos 2 (a - 0) m } m * 2 In ^ + ^(« - P) m ^| 

m 2 H -m A mjj m 2 h -m A m 2 h 

2 m 2 1 

-^om!(— +7E + ln^, (16) 

z - = 1 -(sW{5 (2m ° + " i - + " i ' +m - )+ #i li, 5 

+ cos 2 (a - (3) 9 G H 9 In — £ + sin 2 a - (3) 9 G \ In 
r?% — vtlq m H m L h — m G 

_ 2iVcm?( _l_ + 7£ + ln ^__ 1) }. (17) 

The wave function renormalization constants of the and Z bosons, on the other 
hand, are given by Z z = 1 + 0[M 2 z jv 2 ), and = 1 + 0(M^/v 2 ), and may be set 
equal to unity in our approximation scheme. In summary, the modification factors 



11 



in (13) turn out to be 

(4tt)V 



r <z -, n w ~ i , 1 N c m 2 

'-'mod — i ) °mod — 1 + 7T~TT~5 ' o ■ 



4. The Higgs-Goldstone system 

Being equipped with the machinery of the equivalence theorem, we are now 
interested in the radiative corrections to the processes H — > zz and H — > w + w~. 
Since we assume that all the Higgs boson masses are much greater than M z and 
Mw, internal loops are also dominated by Higgs and Nambu-Goldstone bosons (plus 
top quark). 

Before starting to calculate the radiative corrections, we have to go somewhat in 
detail on the structure of the counterterms, since there are subtleties with regard to 
the field mixing. The interaction Lagrangian relevant to Hzz and Hw + w~ vertices 
is extracted from the Higgs potential (1) 

2 

C Hzz = ^sm(a-P)Hzz, (19) 

2 

C Hww = ™^ sin(a - (3)Hw + nj-. (20) 
v 

The counterterms, 5Chzz, and SChww, consist of two parts, i.e., 

SChzz = SChI z + S^Hzzi 8£hww = 5£hLw + ^Hww ( 21 ) 

The first one is obtained simply by replacing the parameters in (19) and (20) as 
m 2 H — > m 2 H — 5m 2 H , a — > a — 5a, (5 — > (5 — 8(5, v — > v — 5v. We thus obtain 



. , : I 5 1 HI, 5v 1 III], . 

m 



(5a - 5(3) cos(a - (5)Hzz, (22) 
2v 
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5C 



(i) 

Hww 



StTL^ Si) I Tft^ 

' ■ — } —sm(a-/3)Hw + w- 

V V 



m 



^{5a - 5(3) cos(a - (3)Hw + w~ . 



(23) 



The second one is rather complicated due to field mixing effects between H <-> h, 
z ^> A, and w ± <-* G ± pairs. These mixing effects are described by 2 x 2 wave 
function renormalization matrices together with the mixing angle renormalization, 
5a and 5f3. The renormalization is fullfilled by the replacement 

\ / 




cos 5a —sin 5a 
sin 5a cos 5a 

cos 5(3 —sin 5(3 
sin 5/3 cos 5(3 

( cos 5 (3 —sin 5(3 
sin 5/3 cos 5/3 




(24) 



(25) 



(26) 



Because of these mixings, the counterterms for the Hzz vertex are also produced 
from the hzz and HAz terms in the Lagrangian, i.e., 



C-hzz = - ^ cos(a - (3)hzz, 



■-HAz 



m H -m A 



cos(a — (3)HAz. 



(27) 
(28) 



Similarly, the interactions coming from the potential (1) 



-■HGw 



m 



- cos(a — (3)hw + w , 



mH m ° cos(a - (3)(Hw + G- + Hw~G + ), 



(29) 
(30) 



provide some of the counterterms for Hw + w interaction. We will restrict our 
consideration to the lowest-order loop-corrections and hence keep only linear terms 



13 



in 5a and 6(3. It is in fact straightforward by the replacement (24), (25) and (26) 
in the interactions (27)- (30), to reach the following set of counterterms which are to 
be added to (22) and (23); 

8C$„ = {(yfa-l) + (Z,-l)}^8in(a-l3)Hzz 

2 

-{sJz^h - 5a)^co S (a - (3)Hzz 
-(y/zZ + fl?) ™*"^ cos(a - (3)Hzz, 

(31) 

2 

MhL = {(y/^-l) + {Z w -l)}^sm(a-l3)Hw + w- 

2 

-(JzZi - 5a)— cos(a - (3)Hw + nj- 
v v 

v V 

(32) 

It is by now clear what kinds of two-point functions are to be computed in order 
to complete our counterterms. Hereafter we will use the notations Uij(p 2 ) for the 
two point functions corresponding to Fig. 1, where indices i and j denote either one 
of the scalar field, H, h, z, w, A or G. Internal particles in Fig. 1(a) and 1(b) are 
given in Table 2. Calculations of some of the two-point functions are sketched in 
Appendices C and D. The basic quantities in the counterterms are summarized as 
follows: 

5m 2 H = Re(u HH (m 2 H )), Z H = 1 + Re (W HH (m 2 H )) , 

Z z = 1 + Re (14(0)), Z w = l + Re(U' ww (0)). (33) 

Here the renormalization is done on the mass shell, namely, at p 2 = m 2 H for Hhh and 
p 2 = for U zz and H ww . The divergences in the mixing of the two-point function 
UhH are taken care of by the counterterms 5a, \JZ hH and ^/Z Hh . In general, Z hH 
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and Zjjh are not bound to be identical. We will, however, impose ZhH = Zuh just 
for simplicity. If the subtraction is carried out at p 2 = m 2 Hl they are determined by 

5a = ^2 Z —2 Re ( n ^«)) + ^ Re ( n UO) , (34) 

= = l -Re (n; ff K)) . (35) 

In the same way, the remainimg counterterms are determined by 

5(3 = ~^i Re {UzA{0)) ~ ^ Re (n ' 2A(0)) = ~\ Re (n ' 2A(0)) ' (36) 
^Z^=^/zZ=^Re(W zA (0)), JZ~~ G = v /^=^Re(n^ G (0)), (37) 

In Eq. (36), we have used the fact that 11^(0) vanishes as a consequence of the 
Nambu-Goldstone theorem. The renormalization of (3 could be done alternatively 
by using Tl wG (0). The difference between Tl zA (0) and Tl wG (0) is of course finite and 
infinities may be eliminated no matter whichever choice we would take. The use of 
U Z A instead of IT^g to determine 5f3 is just our convention. 

We should add a few words as to the counterterm 5v to the vacuuum expectation 
value. In the on-shell renormalization scheme [23] in which Mz, My/, and the fine 
structure constant (e 2 /47r) are physical input parameters, 5v is determined by 
5v 1 (5M&_ _SMI 5 M 2 - 5M 2 W _ 6f\ 

v 2\M 2 V M 2 ' Ml - M; v #)• [66) 
There is an alternative way which is referred to as modified on-shell scheme [24] and 
uses the muon decay constant instead of M w as an input. This scheme was used 
in our previous work [5]. We have confirmed that the difference between the two 
prescriptions is small numerically. Just for definiteness, we use (38) in our numerical 
calculations. 
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5. Loop corrections to the vertices 



We are now in a position to carry out the loop calculations of the Hzz vertex 
together with Hw + w~'s. The Feynman diagrams contributing to the renormaliza- 
tion of the Hzz vertex are depicted in Fig. 2. They are classified into scalar part 
and top quark part: 

t Hzz (p 2 ) = r£ la V) + N c rt7 k \p 2 )- (39) 

The scalar loop contributions are divided further into three terms 

rE ar) ( P 2 ) = rgL + rgL + rgL (40) 

according to the types of Feynman diagrams in Fig. 2. Internal particles in Fig. 2 
are listed in Table 3. 

The calculation is tedious but straightforward and we just record the results. 
The Feynman diagrams corresponding to Fig. 2(a) are summed up to: 

r£L = -^sin 3 (a- (3)g(p 2 ,0,0,m 2 H ) 

4 2 

mhmH sin(a - 0) cos 2 (a - (3)g(p 2 , 0, 0, m 2 ) 



v 3 



1 f 9 .cos a cos 2 3 sin a sin 3. n , . . 

+ -,{ m 2 H ( , - - 2m 2 A sm «- 0) 

v 3 \ HK sm/3 cos/? ' . 

X \ ( m h ~ m Af sin2 (« - P)g(p 2 , m A, m A, m h) 



+ {m 2 H - m 2 A ) 2 cos 2 (a - 3)g(p 2 , m 2 A , m 2 A , m 2 H ) J 

■ sin(a — 3) cos 2 (a — (3)g(p 2 , 0, m^, m 2 ) 



2m 2 (m 2 H -m 2 A )(m 2 h -m 2 A ) ^ 

,3 



,2 \2 



2m 2 H (m 2 H m 2 A y _ ^ _ ^ 



3m% 9 . _ . cos 3 a sin 3 a . . 9 9 9 

+ ^ sm V_ W (___ )9(p VI,, m? ,,o) 
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2m 2 h m 2 H (2m 2 H + m\) 2 sin 2a 2 2 2 

sin(a - /g) cos 2 (a - /3) , 0/3 g(ft >,m 2 H ,m 2 h ,0) 



v 3 v y v y sin2/^ 

+ 2K-^)K- m i)(2< + m D _ sin2a 
u d sin2p 
xg(p 2 ,m 2 H ,m 2 h ,m 2 A ) 

m ^( m H + 2m D ■ / ,o\ 2/ m sin 2a 2 2 2 

ft/ sin(a - (3) cos 2 (a - (3)- -g(p 2 ,m 2 h ,m 2 h ,0) 



v 3 sin 2/3' 

(ml - m\)\m 2 H + 2m 2 ) ^ _ p) s^2a ^ ^ ^ } 



t> 3 sin 2/?' 

(41) 

Each term of the above expression is easily identified with the corresponding Feyn- 
man diagrams just by looking at the function g(p 2 , m\, m\, m 2 ) defined in Appendix 
B. Those diagrams of the type of Fig. 2(b) are, on the other hand, expressed by the 
integral j 2 (p 2 ,m\,m 2 2 ) which is also defined in Appendix B. They turn out to be: 

IffL = ^jf{mlcos 2 (a- (3) + m 2 H sm 2 (a- /?))sin(a- (3) f 2 (p 2 ,0,0) 



1 f o /Cos a cos 2 /3 sin a sin 2 j3 2 . . . 

m H ( : — — ) — 2m A sm(a — p) 



2v 3 \ ny sin/5 cos (3 

x < 3m\ sin 2 (a — (3) + 3m 2 H cos 2 (a — (3) 



sin 2a , 2 i \ I r / i i i \ 
+ ^^( m h- m H))h(j> ,m A ,m A ) 

2 f 2 cos a cos 2 (3 sin a sin 2 (3 . , 
v 3 \ M sin/? ' cos/5 



) — 2m 2 G sin(a — /3) 



r i i 

x< -m 2 h sm 2 (a — (3) + -m^ cos 2 (a — /3) 

I _ 

1 sin 2a , 2 2 \ , 2 1 /■ / 2 2 2 \ 
2 

~rf( m h ~ m H)( m H - m c) sin (« - P) cos 2 (a - (3)f 2 (p 2 , m G , 0) 
3 

— 3 - m H)( m H ~ m \) sin (" - P) cos 2 (a - (3)f 2 (p 2 , m\, 0) 

3m 2 H \ ( 2 2 \ sm 2« n o\ 2/ m 2] 
' ' (m£ - m 2 H ) , + 2m^ cos 2 (a - /3) + } 



2v 3 \y h " y sin2/3 
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X 



cos 3 a sin 3 a " 



sin (3 cos (3 



hip , m H, m n) 



1 



sin 2a 
sin 2/5 



+ 2m 2 A(m 2 + 2m 2 H ) 



x sin(a - /3) cos 2 (a - /9)^|^/ 2 (p 2 , m^, mj[) 

sin £0 

X I / , o 9 \ Sill 20! 9 \ 9 / 9 

(m h — itl h j— + 2m A sin (a — p) + m ft 



2w 3 



sin 2(3 



(o 2 \ * / zO\ 2Q! . Q o 9 \ 

+ 2m J sm(a - /3)-—-f 2 {p , m h , mj. 
sin 



(42) 



Finally we come to the sum of Feynman diagrams of the type Fig. 2(c): 

1h„ " ^ m ^sin2/5 +2mA j 
x sin(a - (3) cos 2 (a - (3)f 2 (0, m 2 h , 0) 

2mi f. 9 sin 2a „ . , 
+— ^ (m 2 — — + 2m 2 4 ) sm(a - (3) cos 2 (a - (3) 
v 3 [ sin 2p 

, .sin 3 a cos 3 a. 1 „ . 9 

+ro -W-5^ )8n,(<, " /?) } A(<),ro ' ,<)) 

_ 2K- m j,) f g in2a _ 
ir [ sin 2p 

— cos(2a — 2/3) | sin(a — (3)f 2 (0, m 2 A , m 2 ) 

m% — m\ f , 9 sin 2a n 9 . . _ „ 

( m ft~ + 2m 2 A ) sin(2a - 2/5) 



1/ sin 2/3 

2 , cos 3 a sin 3 a 1 m *■ / 2 2 

+ 2r "M^-T o) cos (« - P) f cos (« - P)/2(0, m A , m H ). 

sin p cos p I 



(43) 



The top quark contributions are obtained by using the Yukawa couplings listed in 
Table 1: 

r gua rk) = 4m^cosa ^ ^ _ ^ ^ ^ ^ j (M) 



18 



The loop corrections to the Hw + w vertex go in the same way and we again 
separate the contributions into scalar and top quark parts, 

t Hww (p 2 ) = rEf V) + A^rfcV)- (45 ) 

The scalar contributions ^hww T \p 2 ) were computed in our previous paper (see Eq. 
(31) in Ref. [5]), and we need not reproduce them here. The top and bottom quark 
part is given, if we set m b = 0, by 

rtT } = ^ c ^{h(p 2 ,rn 2 ,rn 2 ) + f 2 (o,o,rn 2 )} 

Arrfi cos a . 9 9 9 . 

-g(p 2 ,m 2 t ,m 2 t ,0). 



v 3 sin (3 ' 

(46) 

By adding the counterterm contributions we end up with the decay width for- 
mulae, 

1 1 AM 2 

T(H - Z l Zl) = - , 1 - ~^\M HZZ (P 2 = ™ 2 H )\ 2 \CU\ (47) 

6lTi ttih \ rriH 



where the invariant amplitudes are given by 

M Hzz (p 2 ) = T Hzz (p 2 ) + ^cos(a-P)Re(u hH (m 2 H )) 
— sin (a — (3) Re (n##(m?r 



AM 2 

]V >M H Up 2 = rn 2 H )\ 2 \CZ d \\ (48) 



v 

-Re (W zA (0) + U' hH (m 2 H )) cos(a - (3), (49) 



M Hw w(p 2 ) = ? Hww (p 2 ) + J cos(a - /3)Re (lL hH (m 2 H j) 
— sin(o; — (3)Re (n##(m^)) 
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w 




cos(a — f3) 



I 




v 



sm(a — P) 



(50) 



These amplitudes are necessarily finite and the finiteness is a non-trivial check of 
the calculations. 

6. Numerical analyses of the decay widths 

Let us now analyze the decay width formulae numerically and look at their heavy 
Higgs boson mass-dependences. In doing so we have to select reasonable numbers 
for the set of parameters, m#, mg, m^, m^, a and (3. There are four kinds of 
experimental information that we have to bear in our mind: (1) the measurement 
of the p-parameter [20, 25], (2) the neutral meson mixings ( K° — K°, B° — B°, and 
D° — D°) [26], (3) the recent measurement of the decays such as B — > K* (892)7 
[27], (4) the ratios, Rb = Y(Z — > bb)/T(Z — > hadrons), and the R c counterpart [28]. 

It has been known rather well that the constraints on the deviation of the p- 
parameter from unity prohibits vtiq to be much larger than or much smaller than 
either of neutral Higgs boson masses. These constraints on the Higgs masses, how- 
ever, depend on the mixing angle, a and (3. The simultaneous analysis on the masses 
and mixing angles would become much involved. For the purpose of getting an in- 
sight into a global picture of the ma- and m^-dependences of the decay widths, we 
may well vary the masses , and rriA, a little beyond the p-parameter constraint. 
The data on neutral meson mixings and b — > 57 decay rate both rule out small value 
of m G and small tan f3. Grant [29] has confronted the two-doublet model with these 
experimental data simultaneously and made an overall analyses. The constraints he 
derived are mc > 150 — 200GeV and tan/3 > 0.7. We will give these values due 
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considerations in the following numerical analyses. The so-called Rb — R c crisis [28] , 
which might jeopardise the standard model, is our recent central concern. It seems, 
however, yet premature to draw definite conclusions from those data or implications 
to the two-doublet model. We will henceforth wait a little while for what comes 
next. 

Besides these sets of experimental information, we have theoretical contraints on 
the masses in the two-doublet model. One of them is the triviality bound [30] and 
another is the tree unitarity bound [31]. The latter constraints lead to m H < 500 
GeV, m G < 870 GeV, m h < 710 GeV and m A < 1200 GeV. The triviality bounds 
also provide roughly similar results. We will take these into our consideration as a 
guide of our parameter choices. We will set m# = 300 GeV, and mh = 400 GeV 
throughout our numerical calculations. 

First of all, we look at the m G - and m^-dependences of T(H — > Z L Z L ) and 
T(H — > W^W£) for the following two cases, 

Case I: tan (3 = 2 and sin 2 (a — (3) = 1, 
Case II: tan/5 = 10 and sin 2 (o; — j3) = 0.5. 

The decay widths T(H — > Z L Z L ) and T(H — > W^W£) as a function of m G depicted 
in Figs. 3 and 4 correspond to the above two cases. We vary m G from 300 GeV to 
1000 GeV, while m A is set equal to 400, 600 and 800 GeV. 

To get a rough idea of numerical values let us recall that these decay rates in the 
minimal standard model (MSM) with a single Higgs doublet give at the tree level 

1 m 3 , 



^L^L) tree 



MSM/tt s w+W-M _ 1 H 



[ (H^W+WE)\ tTee 



327T V 2 \ 

1 m 3 



4M: 



1-^. (51) 



m 



H 



16"7T V 2 \ 

For m H = 300 GeV these formulae give r MSM (F -> Z L Z L ) = 3.5GeV and 
pMSM^ _^ W£W L ) = 7.5GeV. The tree level decay widths in the two Higgs 
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doublet model are suppressed by a factor sin 2 (a — ($) compared with (51). To take 
into account the one-loop corrections in MSM amounts to multiplying (51) by the 
factor 



This factor is about 1.013 for m# = 300GeV and the radiative corrections have been 
said to be very small in MSM. 

The situation in the two-Higgs doublet model could differ very much from this. 
The radiative corrections would be enhanced if m|./t> 2 , m 2 A /v 2 and m 2 h /v 2 are large 
and could stand out from the lowest-order values. Our interest is to what extent 
they could supersede the tree-predictions. Figs. 3 and Fig. 4 give us typical patterns 
of the me- an d m^-dependences. One can immediately notice that the decay rate 
of if — * Z L Z L depends on m G and m A surprisingly little in contrast to our naive 
expectation. We have also confirmed that the rrih dependence is very weak. We can 
hardly see any indication of power-behaviors w.r.t. heavier Higgs boson masses. 

The decay width of H — > W^W£ on the other hand exhibits some power- 
behavior. The one-loop corrections are potentially large to overcome the tree-level 
predictions. One thing that we should pay attention here is that, while the radia- 
tive corrections are large in general, they tend to become small for mc = m A . The 
apparent difference lying between the two decay modes into W- and Z-pairs is puz- 
zling. If the isospin symmetry SU(2)y is exact, these two decay rates differ simply 
by a factor of two (up to the phase space difference) irrespectively of the choice of 
the parameters. In our two-Higgs doublet model, the isospin symmetry is broken 
by the A 5 -term in which we have necessarily to seek for the source of the difference. 

The overall properties of the decay widths may be seen easier if one would draw 
the widths in the birds' eye views. Figs 5 and 6 give those of the decay width of 
F(H — > ZlZl) as a function of mc and m A , and Figs. 7 and 8 those of T(H — > 
W£W£) . Figs. 5 and 7 correspond to the case (I) and Figs. 6 and 8 to the case 




(52) 
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(II). 

Figs. 5 and 6 show that the decay width of H — > Z L Z L does not go far away 
from the tree-level values (3.5 GeV for Fig. 5 and 1.75 GeV for Fig. 6) for a wide 
range of parameters, 400GeV < m G , tua < lOOOGeV. This is rather in contrast with 
our naive expectation. Considering that rn^/v 2 and m 2 A /v 2 are greater than unity, 
the radiative corrections could instabilize the tree predictions. The reason for the 
almost flat behavior in Figs. 5 and 6 will be elucidated qualitatively in Sec. 7. 

In the case of H — > W^W£ , on the other hand, we can immediately see in Fig. 
7 a conspicuous ridge elongated along the line mc = tua- The width comes close 
to the tree level prediction along this line, i.e., the radiative corrections tend to 
be minimized in the custodial SU(2)y symmetric limit. The corrections become 
larger and larger as we go off from this mc = rriA line. This reflects the fact that 
the radiative corrections entail the terms proportional to m G , m 2 G m 2 A and m A that 
exceed the tree-values. The shape of the surface in Fig. 8 corrsponding to the case 
II differs from that in Fig. 7. The one-loop predictions go up and down depending 
on the choice of mass parameters. The point to be emphasized is that the line 
m G = rriA is again given a special meaning. The radiative corrections along this line 
tend to be minimized, the tree-level width being equal to 3.8 GeV. 

To sum up our numerical computation, two questions have emerged naturally. 
One is why the radiative corrections to the decay rate of H — > W^W^ are minimized 
for m G = rriA. The other question is what explains the difference between the two 
decay rates and how. We have included the top quark mass effect in Figs. 3-8. The 
top contributions, however, modify the predictions of the widths only within a few 
per cent and are not very significant. 

7. New screening theorem for Higgs vertex 
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To shed light on the two questions raised at the end of Sec. 6, we disentangle 
the leading terms among M.h ww (p 2 = ™>h) an d M-Hzzip 2 = By leading terms 

we mean those containing m 2 G , m 2 A , and m\, but not m 2 H at all. The results are as 
follows: 

M H ww{.m 2 H ) — > \ sin(a - (3) \ (m 2 G - m 2 A )m 2 A - m 2 G m 2 A In I 

[AtiYv 6 [ m A J 

+ (terms depending on the prescription for 5f3) 
+(9(^,^,or^i), (53) 

M Hzz (<) — or ^). (54) 

The second line in (53) is due to the terms that arise by our choice of H z a instead of 
n„,G in defining the renormalization 5/3, i.e., Eq. (36). In other words, these come 
from the last line of Eq. (50). 

This leading behavior derived by hand in (53) and (54) is of course consistent 
with our numerical computations. The decay into a W-pair behaves as fourth power 
of me and/or ttia, but does not contain nth in (53). It is therefore insensitive to 
mh- Furthermore the leading terms in (53) also disappear if we put mc = mA- 
Recall that the second line in (53) consists of those terms containing m 2 G — m 2 A . The 
behavior (54) indicates that all of the leading terms just disappear after summing 
up all the diagrams together with counterterms and the decay rate T(H — * Z L Z L ) 
is fairly insensitive to ma, ttl a , or m^. 

From the viewpoint of the custodial SU(2) V symmetry and its breaking, it is 
natural to divide the decay amplitudes into two parts 

M Hww = M s + M B Hww , (55) 
Mhzz = M s + M B Hzz . (56) 

Here M s is the custodial symmetric part and M^wm an d -M/Lz are those due 
to the custodial symmetry breaking. As we mentioned in Sec. 2, the custodial 
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symmetry breaking occurs through the As-coupling. This breaking effect shows up 
in two ways in the perturbative calculations. One is in the broken- symmetry in the 
Higgs-Goldstone couplings, the other in the different masses in G and A propagators. 
Both of these two effects are collected in the second terms of (55) and (56). The 
symmetric part M. s are obtained simply by replacing by vtiq everywhere in 
M Hww 

and Mhzz- 

The leading behavior shown in (53) and (54) indicates that the custodial sym- 
metric part M. s does not possess the fourth-power terms w.r.t. mo, tra or mh- We 
may conclude that M s is insensitive to the heavier Higgs boson masses. We are 
thus led to postulate a new screening theorem. 

Theorem: There occurs a cancellation mechanism among the leading terms w.r.t. 
the heavier Higgs boson masses in the custodial S , f/(2)y-symmetric limit in the 
radiative corrections to the decay rates of H — > Z L Z L and H — > W^W^. 

It should be stressed that this theorem differs from Veltman's in that the Velt- 
man's theorem has been proved for the oblique-type radiative corrections, while the 
above one is for the Higgs decay vertex. It should be also noticed that despite this 
difference the custodial symmetry is the key ingredient in both cases. Recall that 
the proof by Einhorn and Wudka relies heavily on the custodial symmetry. What 
we did in Sec. 6 and in (53) and (54) amounts to confirming the above theorem 
explicitly on the one-loop level by numerical calculations and by hand, respectively. 
We do not present here a general proof of our theorem and we should not go so far as 
to say anything definite as to the two-loop and higher-order cases, only mentionimg 
the following: the fact that the custodial symmetry is playing the key role in the 
Veltman's theorem is suggestive of the validity of our theorem beyond one-loop. It 
is also extremely tempting to speculate that there could occur a similar cancellation 
mechanism in other decay modes of the if -boson. We will come to this important 
issue in our future publications. 
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Now let us move on to the custodial symmetry breaking part. The absence of the 
leading terms in A4 S , means that the sensitive properties of the amplitude all come 
from M.jj ww and M^ zz . In the J¥-pair dacay, M.jj ww vanishes for m G = m A that 
corresponds to the flat direction shown in Figs. 7 and 8. Thus the answer to the 
first question raised in Sec. 6, namely, the minimization of the radiative corrections 
along the line mc = mA, becomes almost self-evident once we accept our screening 
theorem. 

The second question is what explains the apparent difference lying between de- 
cays into W- and Z-pairs. As we remarked before, the custodial symmetry breaking 
effects appear in the Higgs couplings on one hand, and the different G and A masses 
in the propagators on the other. Let us concentrate on the former, and in particu- 
lar on the broken-symmetry in the Higgs coupling triggered by the As-term in the 
potential 

A 5 (Inr$]> 2 ) 2 = -^[(w + G~ - G + w~) 

+iA {v + hcos(a — (3) — H $m(a — (3)} 

— iz{hsin(a — (3) + H cos(a — (3)} . (57) 

A close look at this expression shows that there exists an interaction (w + G~ — 
w~G + )A, while another possible term 

{w + G~ - w~G + )z (58) 

is missing. The absence of this triple interaction indicates that there are some 
Feynman diagrams contributing to H — > w + w~ which do not have a conterpart in 
H — > zz. Actually the ^-interactions with charged scalars are peculiar: not only 
(58) but also terms such as 

w + w~z, G + G-z, (w + G- +w-G + )z, (59) 

are absent in the whole of the potential V($i, $2)- These are all forbidden interac- 
tion vertices if either G-parity or CP-invariance is exact. 
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This peculiarity leads to a considerable simplification in tracing the m G depen- 
dence in M.uzz on which we now focus our attentions. Neither of T^ zz nor T^^has 
the m G -dependence at all, as one can make sure of in Eqs. (41) and (43). This 
is simply because there is no Feynman diagram of the type of Fig. 2(a) or Fig. 
2(c) involving G ± . In the heavy Higgs mass limit (to g , m^, m/j ^> van ), the m G - 
dependence of the amplitude of the ZlZ L-decay is governed by the combination 

rgL + \ cos(a - (3)Re{U hH {m 2 H )) - X - sin(a - P)Re(U HH (m 2 H )). (60) 

All the other terms in (49) are sub-leading in the heavy Higgs mass limit . The 
m G -dependence comes about from the Feynman diagrams, Fig. 1(a) and Fig. 2(b) 
in which either of wG, or GG pair is encircling. 

Now notice that the last two terms in (60) join together to produce the two-point 
function connecting the state \H > with the linear combination 

^cos(a- P)\h> --sm(a-P)\H> . (61) 

An important point to be observed here is that the couplings of the combination 
(61) to wG- and GG-pairs are exactly of the same strength with opposite sign as the 
quartic zz{w + G~ +w~G + ) and zzG + G~ couplings, respectively, that contribute to 
Tx\ z . Thus the leading terms containing m G {i.e., m G /t> 3 , m 2 G m 2 A /v z and m G m^/f 3 ) 
all disappear in (60) and thus in M.hzz as well. 

The equality of the strength of quartic and triple couplings as described above 
may be undestood easily if we go over to the so-called Georgi-basis [32], in which 
the combinations 

v + h cos(a — f3) — H sin(o; — (3) + iz, 

h sin(a - (3) + H cos(a - (3) + iA (62) 

are taken from the outset. It is also to be mentioned that the equality of the quartic 
and triple couplings is necessary for the ultraviolet divergences to disappear in the 
combination (60). 
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Once the absence of leading terms of the form m G /v 3 , m G ml/v 3 and m G m\/v 3 
is established in A4^ zz , then it is almost trivial that those proportional to m\/v 3 
and m 2 A m 2 h /v 3 are also absent. This is because the leading terms in M^ zz must 
vanish by setting m G = and this is possible only when those terms are not 
there. 

In this way we are able to disentangle so many Feynman diagrams and to trace 
the origin of the difference lying between M.h W w and M.r Z z- We are now convinced 
how and why the moderate behaviors of T(H — > Z L Z L ) are seen in Figs. 5 and 6. 

8. Summary and Discussions 

In the present paper, we have studied one-loop radiative corrections to the Higgs 
boson (H) decay into longitudinal gauge boson pairs, H — > Z L Z L and H — > W£W£ . 
The two Higgs doublet model has been considered throughout and H is assumed to 
be much heavier than the weak gauge bosons. A particular attention has been paid 
to the non-decoupling effects due to the other Higgs bosons which are assumed to 
be all heavier than H . 

Our numerical analyses on the decay rate H — > W^W£ show that the radiative 
corrections are potentially large. The larger the mass difference m G — m\ becomes, 
the larger the deviation from the tree-predictions turns out to be. The point is that 
the radiative corrections tend to be minimized if m G = m\, for which the custodial 
S77(2)y-symmetry is recovered. This fact has led us to postulate a new screening 
theorem for Higgs vertices with reference to the custodial symmetry, which may be 
regarded as a generalization of the Veltman's theorem. 

Although our generalized screening theorem has been confirmed only at the one- 
loop level for the particular decay processes, this type of theorem would hopefully 
play the role of a working-hypothesis in our future study. It will be of particular 
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interest to see if the custodial symmetry would have any relevance to the screening 
effect of heavy particles in other types of non-decoupling processes. Some examples 
to be examined include triple gauge boson couplings [33] and longitudinal gauge 
boson scatterings, W£W L — > W^W^ and W^W£ — > Z L Z L . It is also important to 
see what happens in models other than the two-Higgs doublet model. We will come 
to these subjects in our future publications. 

The decay rate T(H — > Z L Z L ) has turned out to be unexpectedly insensitive to 
heavier Higgs boson masses. The reason for this was elucidated in detail in Sec. 7. 
It is thus rather difficult to use this decay rate to get any signature of the masses of 
unknown heavy particles. Alternatively, however, by turning the tables around the 
decay H — > ZlZl could be useful to get information on the mixing angle a. Fig. 9 
shows the a-dependence of this decay rate for tan/3 = 2, tan/3 = 5 and tan/3 = 10. 
Measurements of the width with an accuracy on the order of a fraction of 1 GeV 
would enable us to determine the value of a. 

Finally we would like to add a few rather peripheral remarks. As we mentioned in 
Sec. 2, we neglected the term (/i^ 2 $J$ 2 + ^12^2^1) which would break the discrete 
symmetry of the Higgs potential (1). By dropping this term, the heavy Higgs boson 
mass limit is rendered to be the same as the strong quartic couplings (see Eqs. (7)- 
(11)). Thus the non-decoupling effects are expected to be potentially large because 
of the strong couplings. If we would included (/if 2 ^1^2 + A*i2^!^i) as m minimally 
super symmetric models, one may wonder whether the non- decoupling effects could 
be expected or not. In the presence of the mass scale //12 in addition to v, the heavy 
Higgs mass limit does not always imply the strong quartic couplings. There exists 
a limit in which fi 12 , m h , m A and m G are all large while the quartic coulings are 
small. In such a case, the decoupling is expected from the begining and the two 
Higgs doublet model becomes similar to the minimal standard model at low energies 
[34]. 



29 



The non-decoupling effects considered in this paper may be studied by the 
electro-weak chiral Lagrangian approach [35-37]. This method is powerful in the 
minimal standard model for systematic studies of low-energy manifestation of heavy 
particles. Whether this method is useful in the two Higgs doublet model as well is 
yet to be scrutinized and we leave it as an open problem. 
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Appendix A 

The Yukawa couplings of top and bottom quarks are given by 

£ Yukawa = ^-^Ql^Or ^~zQ LlT 2 & J R + h.C, (63) 

v sin p v sin p 

for model I and 

^Yukawa = -^^Ol^ " ^^Ol^^A + h.C., (64) 

v cos p f sin p 

for model II, where Ql = (^l, &l) t - More explicitly, these interactions are expanded 
by putting Eqs. (3)-(6) as follows; 

^Yukawa = -m b bb + C{bbE + C 2 bbh + C^i^bz + Cj)i^^bA 
-m t tt + D\ttH + D 2 ith + D 3 ti-f 5 tz + D 4 ii>y 5 tA 
+E 1 [t{\ + i b )bw + + 6(1 - j 5 )tw-} 
+E 2 {t(l + l5 )bG + + 6(1 - 75)^-} 
+Fi {i(l - 75 )6u>+ + 6(1 + -y 5 )tw~} 

+F 2 {f(l - 75 )6G + + 6(1 + 75)^-} . (65) 
The coefficients in the above are tabulated in Table 1. 
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Appendix B 



We define the following functions to express the various Feynman integrals: 

d D k i 



fi(m ) = [i 



4-D 



(2tt) d k 2 -m 



2 



m 2 ( 2 , m 2 

1 + 7jB + In 



2 ' 



(4vr) 2 \D-4 4vr/i 

(66) 

h(p\mlml) = -^- D J ^ Dk / mUp _ k y_ mr ^) 
gtf, ml ml m 2 ) = J ^ _ _ _ _ _ _ _L _ _ _ 1 _ . (68) 

The function (68) can be expressed in terms of some combinations of the Spence 
function. Details are given in Appendix A of Ref. [5]. 



Appendix C 

Here we would like to summarize the full expressions of the various two-point 
functions used in the text. Let us begin with H zA (p 2 ) and U zz (p 2 ), which are both 
vanishing at p 2 = as a reflection of the Nambu-Goldstone theorem and are put in 
the following forms 

Tl zz (p 2 ) = fl zz (p 2 ) - IMO), Tl zA (p 2 ) = fl zA (p 2 ) - fl zA (0). (69) 

The functions fl zz (p 2 ) and fl zA {p 2 ) are given as a sum of scalar and quark contribu- 
tions 

n zz (p 2 ) = nir lar Hp 2 )+nlr rk) (p 2 ), ti zA ( P 2 ) = niT lar) (P 2 )+nir k) (p 2 )(70) 

and each of these terms is given as follows: 

n(scalar) (p2) = < _ Q) 

V 
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777- 

+ ^fsm 2 (a-[3)f 2 (p 2 ,m 2 H ,0) 



(m h - m A ) 

v 2 

(m 2 H -m 2 A ) 2 



sin (a- (3)f 2 (p ,m h ,m A ) 



■cos 2 (a- (3)f 2 {p 2 ,m 2 H ,m 2 A ), 



ttr rk Xp 2 ) = 2N c ^f 2 ( P 2 ,m 2 ,m 2 ) P 2 , 



(71) 
(72) 



fl ( :f aT \p 2 ) = ml{W \ v2 ^ sin(2a-2/?)/ 2 ( P 2 ,m 2 ,0) 

m 2 H {m\ - m|) ^ _ ^ Q) 



2v 2 

m\ — m 2 A . , _ , 9 . sin 2 /3 cos a cos 2 /5 sin ct . 

+ 2 sin a - (3 {m 2 + 

f 2 cos p sin p 

+2m^ cos(a - (3)}f 2 (p 2 , m 2 h , m 2 A ) 

m 2 N — m 2 A . / _. r 9 . cos 2 /9 cos a sin 2 /3 sin a. 

+ 2 cos a - /3 {m 2 ^ ^— 

v 2 sin p cos p 

-2m 2 4 sin(a - f3)}f 2 (p 2 , m 2 H , m 2 A ), 

(73) 

v 2 

Similarly, U ww (p 2 ) and U w g{p 2 ) are given in the same form as in Eqs. (69) and 
(70). fl^ alar )(p 2 ) and fl^ alar - ) (p 2 ) were give explicitly in Ref. [5] and the top quark 
contributions are 

nl q r k V) = -2N c ^(m 2 -p 2 )f 2 (p 2 ,m 2 ,0), (75) 



nirV) = 2A^cot/3/ 2 (p 2 ,m 2 ,m 2 )p 2 . (74) 



nir k) (p 2 ) = -2iV c ^cot/3(m 2 - P 2 )f 2 (p 2 ,m 2 ,0). (76) 
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Appendix D 



The other two-point functions Uhh(p 2 ) and Ylhnip 2 ) were given in Appendices 
C and D in Ref. [5] except for the quark contributions, 



r4T k) = -2iV c ^^{2/ 1 K) + (4m?-p 2 )/ 2 ( P 2 ,m?,m?)} 

N c ( cos 3 a sin 3 oA (quark) N c sin 2a mT (quaxk) 

(77) 

= -2^^ Sin si a n2 C °^ {2/i (m 2 ) + (4m 2 -j>W,™?,™?)} 

+ A^sin2« j _ ^ + gin(a _ ^ quark)} 
f sin 2p L J 

(78) 

Here the tad pole effects on the H and h fields due to the top quark are denoted by 
^quark) and 7^ quark ) respectively. These are given by 

rr „ k)=4 m?cosa /i(m?)i 
f sin p 

T ( quark ) =4 ^sina 
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Table 1. 



Yukawa couplings in (65). 



Coefficients 


Model I 


Model II 


Ci 


—(mb/v) cos a/ sin/3 


(mb/v) sin a/ cos /3 


c 2 


—(mb/v) sina/ sin /? 


—(mb/v) cosa/ cos /3 


c 3 


-m b /v 


-m b /v 


c 4 


—(rrib/v) cot /3 


(mb/v) tan /3 


D 1 


—(m t /v) cos a/ sin/3 


—(mt/v) cos a/ sin /3 


D 2 


—{m t /v) sina/ sin /3 


—(m t /v) sin a/ sin /? 


D 3 


m t /v 


m t /v 


D 4 


(mt/v) cot /3 


(m t /v) cot /3 


E 1 


-mj / a/2v 


— mfc / a/2v 


E 2 


— (m b / V2v) cot /3 


(mb/\/2v) tan/5 


^1 


m t /\/2v 


m t /V2v 


F 2 


(m t /V2v) cot (3 


(m t /V2v) cot /3 
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Table 2. 



Combinations of internal particles (X, Y) running in Fig. 1 (a) and X in Fig. 
contributing to Ii zz and Ti z A ■ 



Propagator 


type 


Internal particle species 




Fig. 1(a) 


(H,z), (H, A), (h,z), (h,A), (t,t) 




Fig. 1(b) 


G, A, h, H 


n zA 


Fig. 1(a) 


(H,z), (H, A), (h,z), (h,A), (t,t) 




Fig. 1(b) 


G, A, h, H 



Table 3. 

Combinations of internal particles (X, Y; Z) in Fig. 2(a) and (X, Y) in Figs, 
and 2(c) for the vertices T^ zz (i = 1, 2, 3) and r^"^ rk \ 



Vertex 


type 


Internal particle species 


r (i) 

1 Hzz 

r (2) 
1 Hzz 

F (3) 
1 Hzz 

-p (quark) 
1 Hzz 


Fig. 2(a) 

Fig. 2(b) 

Fig. 2(c) 
Fig. 2(a) 


(z,z;h), (z,z;H), (A, A; h), (A, A; H), 
(A,z;H), (A,z;h), (H,H;z), (H, H; A), 
(H,h;z), (H,h;A), (h,h;z), (h,h;A), 
(w,w), (z,z), (G, G), (A, A), (A,z), 
(G,w), (H,H), (H,h), (h,h) 
(z,h), (z,H), (A,h), (A, H), 
(t,t;b), (b,b;t) 
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Figure Captions 



Fig. 1 Three types of Feynman diagrams contributing to the two-point functions 
Uij, where % and j are either of H, h, G or A. Fig. 1 (c) denotes the tad pole 
contributions. Internal particles (X, Y) in (a) and X in (b) are listed in Table 
2. 

Fig. 2 Three types of Feynman diagrams contributing to the Hzz vertex; (a) 

and rfe rk) , (b) r%\ z , and (c) rg^. Internal particles (X,Y;Z) in (a) and 
(X, Y) in (b) and (c) are given in Table 3. 

Fig. 3 The decay width T(H -> Z L Z L ) ((a), (b) and (c)) and T(H -> W£W£) 
((A), (B) and (C)) as a function of m G . The choice of itla is varied as 400 GeV 
((a) and (A)), 600 GeV ((b) and (B)) and 800 GeV ((c) and (C)). Mixing angles 
are determined by tan/3 = 2, and sin 2 (o; — (5) — 1, and the neutral CP-even 
Higgs boson masses are m H = 300 GeV and m h = 400 GeV. 

Fig. 4 The decay width T(H -> Z L Z L ) ((a), (b) and (c)) and T(H -> W^Wf) 
((A), (B) and (C)) as a function of m G . The choice of is varied as 400 
GeV ((a) and (A)), 600 GeV ((b) and (B)) and 800 GeV ((c) and (C)). Mixing 
angles are determined by tan/5 = 10, and sin 2 (o; — j3) = 0.5, and the neutral 
CP-even Higgs boson masses are m# = 300 GeV and m h = 400 GeV. 

Fig. 5 The bird-eye's view of the decay width T(H — > Z L Z L ) as a function of mc 
and rriA- The mixing parameters are determined by tan/5 = 2, and sin 2 (a — 
(3) — 1, and the neutral CP-even Higgs boson masses are ra# = 300 GeV and 
m h = 400 GeV. 

Fig. 6 The bird-eye's view of the decay width T(H — > ZlZl) as a function of 
mo and m\- The mixing parameters are determined by tan/5 = 10, and 
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sin 2 (a — (5) = 0.5, and the neutral CP-even Higgs boson masses are tuh = 300 
GeV and m h = 400 GeV. 

Fig. 7 The bird-eye's view of the decay width T(H — > W^W£) as a function 
of m G and ttia- The mixing parameters are determined by tan/3 = 2, and 
sin 2 (a — P) — 1, and the neutral CP-even Higgs boson masses are m H = 300 
GeV and m h = 400 GeV. 

Fig. 8 The bird-eye's view of the decay width T(H — > W^W£) as a function 
of mo and m^. The mixing parameters are determined by tan/3 = 10, and 
sin 2 (a — (3) = 0.5, and the neutral CP-even Higgs boson masses are m H = 300 
GeV and m h = 400 GeV. 

Fig. 9 The mixing angle a dependence of the decay width T(H — > Z^Zl). The 
mixing angle f3 is determined by (a) tan/3 = 2, (b) tan/3 = 5 and (c) tan/3 = 
10, respectively. The Higgs boson masses are m H = 300 GeV, m h = 400 GeV, 
m G = 500 GeV and m A = 600 GeV. 
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